Self-assembly of colloidal polymers via depletion-mediated lock and key binding 
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We study the depletion-induced self-assembly of indented colloids. Using state-of-the-art Monte 
Carlo simulation techniques that treat the depletant particles explicitly, we demonstrate that colloids 
assemble by a lock-and-key mechanism, leading to colloidal polymerization. The morphology of the 
chains that are formed depends sensitively on the size of the colloidal indentation, with smaller values 
additionally permitting chain branching. In contrast to the case of spheres with attractive patches, 
Wertheim's thermodynamic perturbation theory fails to provide a fully quantitative description of 
the polymerization transition. We trace this failure to a neglect of packing effects and we introduce 
a modified theory that accounts better for the shape of the colloids, yielding improved agreement 
with simulation. 



The goal of self-assembly is to tailor the interac- 
tions among nano-scale particles so that they sponta- 
neously assemble themselves into functional materials or 
devices [l] [2]. Such processes are widespread in biol- 
ogy, where they have been optimised by evolution so 
that assembly is rapid and reliable. However, mimick- 
ing this behaviour in the laboratory involves many chal- 
lenges, particularly the design and synthesis of particles 
whose interactions can be accurately predicted and con- 
trolled. Notable experimental successes have included 
assembly of unusual crystals from either "patchy" or 
DNA-functionalised colloids [3] [4] . More recently, par- 
ticles have been shown to self-assemble into structures 
that depend strongly on their geometrical shapes [M9], 
and the role of shape and packing effects in self-assembly 
has also attracted theoretical and computational inter- 
est [7HT2], Here, we use computer simulations to show 
how self-assembly of indented colloidal particles can be 
tuned by subtly varying their shape and interactions, in 
a manner that should be accessible in experiment [6 . 

To this end, we exploit depletion forces [13], which en- 
able the precise control of particle interactions that is re- 
quired for self-assembly. Depletion is an attractive inter- 
action between colloid particles that arises when they are 
mixed with much smaller 'depletant' particles, for exam- 
ple polymers or another species of colloid. These forces 
are particularly strong in colloids with complementary 
geometrical forms, such as buckled spheres [6 or bowl 
shapes [7 . Such systems can assemble via "lock and key 
binding" in which the convex part of one particle inter- 
locks with the concave part of another [6j [TTJ [14] . Fig. IT] 
shows the results of a computer simulation, where col- 
loids with self-complementary shapes [6] have assembled 
themselves into chains, in the presence of depletant parti- 
cles. We show in the following that the properties of these 
"colloidal polymers" can be controlled through the col- 
loidal shape and the depletant number density. The per- 
sistence length of the polymer depends on the colloidal 
shape, and for some shapes, the chains can also branch, 
leading to the possibility of "empty liquid" states [15] [16]. 

We emphasize that the complementary shapes of col- 
loidal particles [6] [17] [18] and properties of the depletion 



FIG. 1: (Color online) Equilibrium snapshot of self- assembled 
chains of spherically-indented colloids. The depletant is 
shown explicitly in the left panel but suppressed in the right 
panel, for clarity. System parameters are h — 0.7 and 
rf s = 0.094: see text for details. 



interaction can both be measured and controlled in ex- 
periments. Indeed, some depletants even allow colloid in- 
teractions to be tuned in situ [5, 19, 20 , potentially lead- 
ing to real-time adaptive control of interaction parame- 
ters [21]. However, the experimental parameter space as- 
sociated with mixtures of colloid and depletant particles 
is very large, depending on the size, shape and concen- 
tration of both species. Theory and computer simulation 
can therefore offer guidance for experiment, by predict- 
ing the parameters for which robust assembly occurs, and 
the likely nature of the self- assembled products. We ar- 
gue that such simulations should deal explicitly with the 
depletant particles, both in the interests of reproducing 
the experimental reality and for avoiding the need to de- 
velop effective ("depletion") potentials, which for irregu- 
larly shaped particles represents a formidable task. 

Thus, we have deployed state-of-the-art Monte 
Carlo (MC) simulation techniques to study spherically- 
indented colloids, together with smaller hard sphere par- 
ticles which act as a depletant. The shape of each in- 
dented colloid begins as a hard sphere of diameter a\ — 1, 
from which an indentation is formed by cutting away a 
sphere of the same diameter, whose center is a distance 
d c from the center of the original sphere. Thus, the di- 




FIG. 2: (Color online) GCA move set for indented particles, 
(a) A colloid reflects through a point pivot to a new position 
(outlined). Any particles that overlap in the new position re- 
flect through the same pivot to occupy the space vacated by 
the colloid, (b) Constrained plane reflections (see text) allow 
for small scale vibrations of colloids within a chain, (c) Re- 
flecting a colloid in a plane passing through its center moves 
it to a new position (outlined). When overlapping particles 
are similarly reflected, the chain 'flexes'. 



mensionless depth of the indentation is h = 1 — d c /ai. 
Our systems contain N = 60 colloid particles in a box 
of size V , with a number density p = N/V = 0.2cr^ 3 
and we consider values of h between 0.3 and 0.7. The 
hard spheres comprising the depletant fluid have diam- 
eter <j s = O.lcji. The colloid shape and the size ratio 
between colloids and depletant are consistent with ex- 
perimental studies such as |B]. 

To obtain accurate computational results for this sys- 
tem, we use a variant of the geometrical cluster algorithm 
(GCA) [22]. This is a sophisticated Monte Carlo scheme 
that samples the Boltzmann distribution of a system by 
updating large groups ("clusters") of particles, with both 
colloids and depletants moving together. Use of such a 
specialized technique is essential for coping with the dis- 
parity in size between colloids and depletant. Standard 
Monte Carlo and Molecular Dynamics techniques are un- 
equal to the task of relaxing such systems because the de- 
pletant acts to frustrate colloidal motion except on very 
small length scales. This problem can be readily appre- 
ciated from Fig. [I] 

The GCA is based on self-inverse geometric operations 
that can be tailored to effectively sample the system of in- 
terest. In the case of self- assembled structures it is essen- 
tial that relaxation occurs on all length scales to ensure 
ergodicity. To this end we use the combination of updates 
described in Fig. ^1 A "pivot" (point reflection) opera- 
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FIG. 3: (Color online) (a) Simulation estimates of X(rf s ) for 
h = 0.3, 0.5, 0.7; lines are guides to the eye and uncertain- 
ties are comparable to the symbol sizes. Comparison of the 
simulation data (Sim) with the predictions of TPT and SCW 
(see text) in the non-branching regime (b) h = 0.7, and (c) 
h — 0.5. There are numerical uncertainties in the TPT/SCW 
predictions which are comparable with those in the simulation 
data: these arise from the numerical estimation of /a from 
simulations containing two colloids with depletant. 



tion (see Fig.l^a)) is employed to relax particle positions, 
while a plane reflection operation [23] allows colloids to 
sample different orientations. For the latter operation, 
if the reflection plane is constrained to lie close to the 
orientation vector of a monomer in the chain, then this 
monomer samples small variations in position and ori- 
entation within its binding pocket (Fig. [2Tb)). On the 
other hand, if the plane is placed through the center of 
the monomer at an arbitrary angle to the orientation 
vector, then large scale flexing of the polymer occurs 
(Fig. [2fc)). Pivot moves in the GCA are reject ion- free 
by construction. However, to ensure the local charac- 
ter of the updates shown in Fig. J2^b), we reject moves 
in which the colloid in its new position overlaps with a 
second colloid. For the move of Fig. |2|c) we ensure con- 
sistency with periodic boundary conditions by rejecting 
the proposed update if the linear dimension of the cluster 
exceeds half the box size. All updates exploit a highly 
efficient hierarchical overlap search algorithm that allows 
us to determine whether a proposed move leads to over- 
laps of our anisotropic particles [24-26 . 

The depletant particles are treated grand-canonically 
in our simulations. That is, their number is free to fluctu- 
ate, corresponding to the common experimental situation 
of a depletant that is in equilibrium with a bulk reservoir. 
We therefore quote the depletant reservoir volume frac- 
tion rf s = N s tt(7 3 /(6V) as a measure of the driving force 
for depletion induced assembly, where N s is the average 
number of depletant particles. 

Turning to our simulation results, we first assess how 



the degree of polymerization depends on the depletant 
volume fraction rf s . To this end, we label the indentation 
on each colloidal particle as its "lock site" . The concave 
part of the surface acts as the "key" , which fits snugly in 
the lock. Let N\, be the average number of lock sites that 
are available for binding (where no other colloidal parti- 
cle is already bound), and let N-& be the average number 
of colloidal particles that are not currently occupying any 
lock site (occupation of a lock site is decided on the ba- 
sis of a radial cutoff criterion; results are insensitive to 
the choice of this cutoff) . The number densities of such 
particles are then pl = N^/V and px = Nk/V. We also 
define X = pi,/ p. For an unassociated fluid, X « 1; for a 
system consisting of long colloidal polymers then X « 0. 
If the polymers are 'tree-like' (without closed loops), then 
the average degree of polymerization is 1/X. 

Fig. u3[a) depicts our measurements of X (black cir- 
cles) for various h, as the depletant volume fraction rf s 
is increased. The range of rj r s over which polymeriza- 
tion occurs is quite narrow in each case (particularly for 
deep indentations), and this range is shifted to smaller 
r] r s as h increases. Physically, the lock- and- key binding 
is strongest when the colloid indentations are deep, and 
the shape complementarity is most pronounced. At the 
largest values of r] s , almost all the colloids are members 
of chains, X w 0. 

Fig.|4]shows snapshots of the equilibrated polymer con- 
figurations that form for h = 0.5 and h = 0.3, at values of 
rf s corresponding to X « 0.1. Compared with the results 
for h = 0.7 (Fig.JT]) one observes that deeper indentations 
result in stiffer chain conformations. To quantify these 
differences, we have measured the persistence length 6, 
defined through (cos(Ol)) = e -L / 6 , where 0l is the an- 
gle between orientation vectors of colloid particles that 
are Lth neighbours in the chain. Thus, large values of b 
correspond to stiff chains: for h = (0.3,0.5,0.7), we find 
b = (1.0, 3.3, 9.1) a \. This wide range of b illustrates how 
the properties of self- assembled colloidal polymers may 
be controlled through the colloid shape. 

For shallower indentations the assembled polymers 
may support branching. This is only possible when the 
indentation is small enough for two colloid particles to 
"lock onto" the convex part of a third one. The marginal 
case is h = 0.5, for which a single key surface can just 
accommodate two locks. As h decreases, the branch- 
ing probability increases rapidly. When bonds are strong 
(X < 0.1 as in Fig. [3, we find that the fractions of 
particles involved in branching for h = (0.5,0.4,0.3) are 
(1%, 9%, 15%). Again, by changing the colloid shape, the 
self-assembled chains can be varied from linear polymers 
(h = 0.7) to chains with a few branches (h = 0.5), and 
finally (h = 0.3) to strongly branched structures reminis- 
cent of the network fluid (or "empty liquid" ) states found 
in studies of "patchy" colloids p~5j [16] . 

Given that this range of behaviour is possible on vary- 
ing just the depletant density rf s and the indentation 




FIG. 4: (Color online) Snapshots for X < 0.1 (depletant not 
shown) and varying indentation depth, (a) At h = 0.5, 
r]l = 0.105, this system consists of just two large chains. 
A single branch point is also indicated, (b) At h = 0.3, 
rf s — 0.14, the polymers form an interconnected network of 
chains, reminiscent of an empty liquid |15| . 



depth /i, theoretical insight is very valuable in guiding 
choices of colloidal geometry and depletant parameters, 
both in simulation and, potentially, in experiment. We 
have applied Wertheim's theory of associating fluids [27] 
to these indented colloids, following the work of Sciortino 
and co-workers [15, 28 on 'patchy' colloids. This theory 
generalizes liquid state theory, incorporating steric con- 
straints. For example, at most one particle may occupy 
any lock site; we also assume that chain branching may 
not occur, which is valid for h > 0.5. Within the the- 
ory, depletion interactions appear as two-body effective 
interactions between the colloidal particles, obtained for- 
mally by integrating out the depletant fluid. Based on 
these assumptions, Wertheim's theory gives a diagram- 
matic series for the density functional of the system, from 
which the number densities p\, and px may be derived. 

At leading order, the theory reduces to the familiar law 
of mass action: plk = PlPkKo where plk = P~ Pl is the 
number density of bonds (i.e. the number density of occu- 
pied lock sites), and Kq is the bare equilibrium constant, 
which depends on the attractive forces between parti- 
cles. The law of mass action applies quite accurately in 
the dilute limit paf <C 1, but to go beyond this limit, one 
must also take account of repulsive forces between parti- 
cles, and the resulting packing effects. The second part 
of Wertheim's theory achieves this by a perturbative ex- 
pansion about a reference system without any attractive 
interactions (rf s = 0). The theory is therefore accurate 
if the packing of particles in the presence of attractive 
forces is very similar to their packing in the reference sys- 
tem. Formally, the thermodynamic perturbation theory 
(TPT) of Wertheim approximates the density functional 
of the system by an infinite subset of terms in its diagram- 
matic expansion. The result is that the bare equilibrium 
constant Kq in the law of mass action is replaced by [28] 

K = n dri2 da;2 # R ( r '12^1^2) /a (^ '12,^1,^2)- (1) 



Here W\,Wi represent the orientations of two particles, 
with ri2 the vector between them, and Q is the phase 
space volume associated with one particle's orientation. 
Also, #R,(r,u;i,u;2) is the two-particle distribution func- 
tion in the reference system (without attractions), while 
fx(r 1 u)i 1 u)2) is a Mayer-/ function associated with the 
attractive part of the effective interactions between par- 
ticles. In the dilute limit, #r = 1 except when the two 
particles overlap, and one recovers the standard formula 
for the bare equilibrium constant K = Kq. Outside the 
dilute regime, particle packing effects are taken account 
of through #r. 

For the spherical patchy particles considered by 
Sciortino and co-workers [28 , #r can be approximated 
from Percus-Yevick theory, and /a is known exactly. 
Thus, the TPT calculation can be performed analytically, 
it describes their data very accurately. For indented col- 
loids, neither #r nor /a is known exactly, but #r may be 
obtained from a simulation of the reference system of in- 
dented colloids in the absence of depletant, and f\ from 
a simulation of two particles in the presence of depletant. 
The integral in (IT]) can then be calculated. The result- 
ing TPT predictions are shown in Fig. [3^b) and (c) for 
h = 0.7 and h = 0.5 respectively: the agreement is rea- 
sonable but there are significant deviations. We attribute 
these deviations to differences in the packing properties 
of chains of particles, compared to isolated monomers. 
(We have ruled out the alternative possibility that many 
body colloidal interactions are responsible by simulating 
indented colloids with prescribed two-body attractive in- 
teractions, for which similar deviations from TPT are 
found [29].) To explore packing effects, we return to the 
diagrammatic analysis of Wertheim, but instead of fol- 
lowing the TPT, we consider just a few simple terms in 
the density functional: see Fig. [5] Under this approxi- 
mation, the law of mass action is replaced by 



Plk = PlPkK [1 + pvi + pLKV 2 ]- 



(2) 



Here, v\ and v 2 are geometrical factors (independent of 
r] r s ) that account for packing of free particles and short 
chains: the relevant liquid-state diagrams are shown in 
Fig. [5] Within Wertheim's TPT, the ^i-term is included, 
but the ^2-term is absent. Further, comparison between 
pi and the law of mass action shows that the effective 
equilibrium constant Plk/(plPk) now depends on the de- 
gree of polymerization of the system, via the ^-term. 
This constant (and hence the quantity X) must there- 
fore be determined self-consistently by solving (J2|, so we 
refer to the analysis in the presence of the v 2 -teiva as 
self-consistent Wertheim (SCW) theory. 

We have obtained values for v\ and v 2 by simulating 
very small systems (up to 4 indented particles, with- 
out depletant). Results are shown in Fig. [5|e). Using 
these values in pi) leads to the SCW predictions shown 
in Fig. [3] In terms of the density functional, SCW theory 
is a much cruder approximation than TPT. However, the 
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FIG. 5: Liquid state theory diagrams showing contributions 
to the density functional that are relevant for calculating /9lk 
within Wertheim theory. Directed heavy lines correspond to 
attractive lock-and-key binding while thin lines correspond 
to repulsive interactions. Where diagrams include dashed 
lines, this indicates a sum over diagrams both with and with- 
out these repulsive interactions. Vertices have weights corre- 
sponding to various combinations of p, />l,Pk, as prescribed 
by Wertheim's theory, (a) Diagram for Ko\ (b) Diagram for 
the vi-term in pi); (c) A diagram included in the TPT but not 
in SCW; (d) Diagrams for the t^-term in Eq. p|, included in 
SCW theory but not in TPT. The terms in (a,b,d) are those 
considered in the SCW theory, (e) Table showing values of v± 
and V2 for various h, in units having 07 = 1. 



SCW theory performs significantly better than TPT, ac- 
counting for around half of the deviation between theory 
and simulation. 

The origin of this effect is the i>2-term in pi). Physi- 
cally, the ^i-term in that equation (included in both TPT 
and SCW theories) reflects the increased virial pressure 
in the system as the colloid number density increases, 
and enhances polymerization. The ^2-term reflects differ- 
ences in packing properties between free particles and as- 
sembled chains. We find v 2 < 0, indicating that as poly- 
merization occurs, the virial pressure is reduced (com- 
pared with TPT), suppressing further chain formation. 
For spherical patchy particles (h = 0), Fig. p)Je) shows 
that the magnitude of v 2 is small, consistent with the 
success of TPT in that case [28]. However, the size of v 2 
grows as h is increased, leading to the deviations from 
TPT shown in Fig. |3| The SCW theory accounts for 
some of these deviations, although agreement is still not 
perfect. One could improve the theory by including fur- 
ther terms in the approximate density functional, but our 
main point here is that packing effects of non-spherical 
particles can significantly affect self-assembly, and are re- 
quired for quanitative theoretical analysis. 

To summarize, the sophisticated Monte Carlo tech- 
niques that we have applied here show that indented col- 
loids can assemble into chains, with persistence lengths 
and branching properties controlled by a combination 
of colloid shape and depletant density. Two variants of 
Wertheim's theory have been analysed, showing how par- 
ticle shapes can influence their self-assembly. Given that 
depletant parameters and the shapes of indented colloids 
can both be controlled in experiments [5j [6j [T7H20] , we 
hope that these results will stimulate further experimen- 
tal studies of self-assembly in these systems. 
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